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Abstract 



This work is devoted to the decay of random solutions of the unforced 
Burgers equation in one dimension in the limit of vanishing viscosity. 
The initial velocity is homogeneous and Gaussian with a spectrum 
proportional to A;" at small wavenumbers k and falling off quickly at 
large wavenumbers. 

In physical space, at sufficiently large distances, there is an "outer 
region", where the velocity correlation function preserves exactly its 
initial form (a power law) when n is not an even integer. 

When 1 < n < 2 the spectrum, at long times, has three scaling re- 
gions : first, a region at very small ks with a time-independent con- 
stant, stemming from this outer region, in which the initial conditions 
are essentially frozen; second, a k'^ region at intermediate wavenum- 
bers, related to a self-similar ly evolving "inner region" in physical 
space and, finally, the usual region, associated to the shocks. The 
switching from the |A;|" to the k"^ region occurs around a wavenumber 
ks{t) oc t~i/[2(2-"')]^ while the switching from k"^ to occurs around 
kiit) oc (ignoring logarithmic corrections in both instances). 

When — 1 < n < 1 there is no inner region and the long-time 
evolution of the spectrum is self-similar. 

When n = 2, 4, 6, . . . the outer region disappears altogether and the 
long-time evolution is again self-similar. For other values of n > 2, 
the outer region gives only subdominant contributions to the small-A; 
spectrum and the leading-order long-time evolution is also self-similar. 

The key element in the derivation of the results is an extension of 
the Kida (1979) log-corrected 1/t law for the energy decay when n = 2 
to the case of arbitrary integer or non-integer n > 1. A systematic 
derivation is given in which both the leading term and estimates of 
higher order corrections can be obtained. It makes use of the Balian- 
Schaeffcr (1989) formula for the distribution in space of correlated 
particles, which gives the probability of having a given domain free of 
particles in terms of a cumulant expansion. The particles are, here, the 
intersections of the initial velocity potential with suitable parabolas. 
The leading term is the Poisson approximation used in previous work, 
which ignores correlations. 

High-resolution numerical simulations are presented which support 
our findings. 
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1 Introduction 



The Burgers equation 

dtv + vdxV = vci^v, (1) 

describes a variety of nonlinear wave phenomena arising in the theory of 
wave propagation, acoustics, plasma physics and so on (see, e.g., Whitham 
1974; Gurbatov, Malakhov & Saichev 1991). It was originally introduced by 
Jan M. Burgers (1939) as a toy model for turbulence. On the one hand, it 
shares a number of properties with the Navier-Stokes equation : the same 
type of nonlinearity, of invariance groups and of energy-dissipation relation, 
the existence of a multidimensional version, etc. On the other hand, it is 
known to be integrable (Hopf 1950; Cole 1951) and therefore lacks the prop- 
erty of sensitive dependence on the initial conditions (chaos). When this was 
realized, interest in the Burgers equation as a toy model for Navier-Stokes 
turbulence abated somewhat. More recently, it became clear that a number 
of questions connected to the Burgers equation are far from trivial. One class 
of problems, arising e.g. in surface growth (Barabasi & Stanley 1995), leads 
to a Burgers equation with a random driving force (usually white-noise in 
time) added to the right hand side (Kardar, Parisi & Zhang 1986; Bouchaud, 
Mezard & Parisi 1995; Cheklov & Yakhot 1995; Polyakov 1995; Gurarie 
& Migdal 1996; E et al. 1997). Even in the absence of a driving force, 
the explicit solution provided by the Hopf-Cole method leads to non-trivial 
problems in the limit of vanishing viscosity, when random initial conditions 
are assumed. An example, motivated by work on large-scale structures in 
the Universe (Zeldovich 1970; Gurbatov & Saichev 1984; Gurbatov, Saichev 
& Shandarin 1989; Shandarin & Zeldovich 1989), is when the initial velocity 
field is a Brownian (or fractional Brownian) process in the space variable 
(She, AureU & Frisch 1992; Sinai 1992; Vergassola et al. 1994; Aurell et al. 
1997). Other examples arise in the sudy of disordered systems (Bouchaud & 
Mezard 1997). 

Another instance, which is at the core of the present study, is the long- 
time behavior of freely decaying solutions in the limit of zero viscosity and, 
in particular the law of decay of the energy. Before reviewing this issue it 
is useful to put it into a broader perspective, namely the law of decay for 
three-dimensional Navier-Stokes turbulence. 

Karman & Howarth (1938; see also Lin & Reid 1963), investigating the 
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decay of high-Reynolds number, homogeneous isotropic three-dimensional 
turbulence, proposed a self-preservation (self-similarity) ansatz for the spa- 
tial correlation function of the velocity : the correlation function keeps a fixed 
functional shape; the integral scale L(t) and the root-mean-square velocity 
u{t) have a long-time power-law behavior; there are two exponents which can 
be related by the condition that the energy dissipation per unit mass should 
be proportional to / L. But an additional relation is needed to actually de- 
termine the exponents. Kolmogorov (1941) realized this and proposed using 
the so-called "Loitsyansky (1939) invariant" as the additional relation and, 
thereby, derived a law of energy decay v?'{t) oc t^^^/^. Proudman & Reid 
(1954) and Batchelor & Proudman (1956) showed that the integral consid- 
ered by Loitsyansky is actually not invariant. Under the assumption that 
the energy spectrum E{k,t) oc k'^ at small wavenumbers k, the Loitsyansky 
integral is, within a numerical factor, just the coefficient of A;*^. For isotropic 
turbulence, the contribution of the nonlinear term to the time derivative of 
the energy spectrum, which is called the transfer function, should be propor- 
tional to fc^ at small k, thereby preventing the constancy of the Loitsyansky 
integral. However, if initially E{k) oc |/c|"' with a "spectral exponent" n < 4, 
then the coefficient of \k\^ will be invariant. This is called here the princi- 
ple of "permanence of large eddies" (PLE). Kolmogorov's argument is easily 
adapted to that case and gives then M^(t) oc jf;-2(n+i)/(n+3) (^^qq^ g_g_^ Frisch 
1995, Section 7.7). 

Almost exactly the same arguments may be applied to the one- dimensional 
Burgers turbulence. In one dimension the transfer function should be pro- 
portional to k"^ at small k (see Section H). Hence, if the initial energy spec- 
trum Eo{k) oc |/c|" with — 1 < n < 2, Kolmogorov's argument predicts again 
u^{t) oc i-2(n^+i)/("-+3)_ shall see that this is correct only when — 1 < < 1. 

For Navier-Stokes turbulence, a rigorous derivation of the law of energy 
decay seems ruled out at present .111 The situation is of course more favorable 
for Burgers-type turbulence and it will be our main goal in this paper to 
critically examine the issues of self-similarity and energy decay at long times. 
In this paper we shall exclusively consider the behavior of the solutions in the 
limit of vanishing viscosity for Gaussian initial data. For the effect of a finite 

^It is not even clear, in three dimensions, if for the initial value problem there is any 
energy decay at all in the limit of vanishing viscosity, since the issue of breakdown of 
smoothness for the Euler flow is moot (see, e.g., Frisch 1995, Section 7.8). 



4 



(small) viscosity and/or non-Gaussian initial data (e.g., piecewise constant 
or shot noise velocities), the reader is referred to recent papers by Funaki, 
Surgailis & Woyczynski (1995), Surgaihs (1996, 1997) and Hu & Woyczynski 
(1996). 

The paper is organized as follows. In Section |^, we formulate our problem 
and list some elementary results about the Burgers equation to be used in 
what follows. Section ^ is devoted to phenomenological derivations of the 
main results about the decay of random solutions to the Burgers equation, 
the nature of which depends crucially on the spectral exponent n. The most 
delicate results, obtained when n > 1, are rederived by systematic methods 
in Sections ^ and |[ In Section H we discuss previous work on this issue. 
In Section ^ we derive the law of decay of energy as the leading term in 
a systematic expansion. Our basic tool is the cumulant expansion (Balian 
& Schaeffer, 1989) of the probability of finding regions in space without 
particles. In Section ^ we calculate the spatial correlation function, and 
prove the permanence of large eddies. In Section |^ we summarize our main 
findings and derive consequences, such as the absence of globally self-similar 
decay for the energy spectrum when 1 < n < 2. Section ^ presents high- 
resolution numerical experiments which confirm such findings. Section ^ 
presents concluding remarks. The Appendix contains a rederivation of the 
Balian-Schaeffer formula in a form appropriate for our problem. 

2 Formulation and elementary results 

We shall be concerned with the initial- value problem for the unforced Burgers 
equation in one dimension : 



Use will also be made of the conservation form of the Burgers equation : 





(2) 




(3) 



and of the formulation in terms of the (velocity) potential : 



(4) 
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where 

V = -d^ijj. (5) 

As is well-known, the solution to the Burgers equation with positive viscosity 
ly, has an explicit integral representation obtained by Hopf (1950) and Cole 
(1951). We are mainly interested here in the solutions in the limit of vanish- 
ing viscosity. Use of Laplace's method then leads to the following "maximum 
representation" for the potential (Hopf 1950) in the limit of vanishing vis- 
cosity (henceforth always understood) : 

'ip{x,t) = max$(a;, a, t) (6) 

(x — a)^ 



$(x,a,t) = 



^o(a) - 



2t 



(7) 



where ipo{0') is the initial potential. Expansion of the square shows that (|^) 
is basically a Legendre transformation. From (|) and (||) it follows that 

v{x,t) = ^^^^, (8) 

where a{x,t) is the coordinate at which $(x, a,t) achieves its (global) maxi- 
mum for given x and t. It is easily checked that a is the Lagrangian coordinate 
from which emanates the fluid particle which will be at x at time t (see, e.g., 
Gurbatov et al. 1991). 

In what follows we shall work with homogeneous random velocity fields. 
We list the main definitions. The correlation function is 

B,{x,t) = {v{x,t)v{Q,t)), (9) 

where the angular brackets denote ensemble averages. Its Fourier transform, 
the energy spectrum is 

1 r°° 

E{k,t) = — B^{x,t)exp{ikx)dx. (10) 

The initial values of the correlation function and of the energy spectrum are 
denoted BQy{x,t) and E^ik), respectively. The (mean) energy is 

/oo 
E{k,t)dk. (11) 
-oo 
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The initial energy is denoted 

/oo 
Eo{k)dk. (12) 
-oo 

We now introduce the basic assumptions of this paper. The initial velocity 
field is homogeneous and Gaussian with an energy spectrum of the form 

Eo{k) = a2|fc|"6o(^), (13) 

where n is the spectral exponent. The even and non-negative function bo{k) 
is assumed to satisfy 60 (0) = 1. It is smooth and decreases faster than any 
inverse power of k at infinity (this implies, of course, that the initial veloc- 
ity correlation function is very smooth). 60 (^) can be characterized by a 
wavenumber ko around which lies most of the initial energy and which is, 
in order of magnitude, the inverse of the initial integral scale Lq. The usual 
definition of the integral scale, which involves the integral of the correla- 
tion function, is not appropriate since this integral vanishes for any n > 0. 
Instead, we use 

Lo ^ ^, (14) 
where is the variance of the initial potential 

ai^mx))= j^^^dk. (15) 

This definition of Lq will be convenient in deriving various long-time laws in 
Sections ^ and |^. 

As we shall see, the value of the spectral exponent n determines the 
universality class of the decaying solution at long times. In order for the 
velocity to be homogeneous its variance must be finite; hence, we must have 
n > —1. When —3 < n < — 1 the velocity has infinite variance and only 
increments are homogeneous; the velocity itself behaves, at large distances, as 
a Brownian or fractional Brownian motion process, a case considered by She 
et al. (1992), Sinai (1992), Vergassola et al. (1994), which is not within the 
scope of the present paper. When —1 < n < 1 the velocity is homogeneous 
whereas the potential behaves, at large distances, as a Brownian or fractional 
Brownian motion process. This case has already been investigated by Burgers 
(1974), by She et al. (1992) and by Avellaneda & E (1995) and will be 
considered here only briefiy. We shall concentrate on the case n > 1, when 
both the velocity and the potential are homogeneous. 



7 



2.1 Elementary results 

We observe that the smoothness of the energy spectrum at the origin depends 
on the spectral exponent n. When n is an even integer, our assumption 
imphes that the initial correlation function BQy{x), being the nth derivative of 
a smooth function of x decreasing rapidly at infinity, also has these properties. 
When n is odd or non-integer the spectrum is not smooth. It follows that 
the correlation function decreases as a power law (Pitman 1968; Bleistein & 
Handelsman 1975), mzi 

Bq^{x) ~ a^Cnlxl""'"^ |x| oo, (16) 

where 

TTTl 

Cn = -2T{n + l) sin—. (17) 

Next, we observe that, when n > 1, the initial potential ipo{x) is a homoge- 
neous random function while, for n < 1, it only has homogeneous increments. 
In the former case, the potential has a translation-invariant correlation func- 
tion, denoted by Bq^{x), the Fourier transform of which is Eo{k)/k'^. In the 
latter case, it is simpler to characterize the initial potential by its structure 
function 

s,^{x) = {iMx)-Mo)f), (18) 

which is still translation-invariant and exists as long as n > —1. It is related 
to the energy spectrum by 



oo 



So4x) = Aa' / (1 - cos(A;x)) P"^ 6o(A;) dk. (19) 







When n > 1 the structure function goes to a finite limit at large separations 
while, when n < 1, it grows without bound: 



Sq^{x) ~ <^ ^2uii-t ^ / 1 |x| ^ oo. (20) 




An important consequence of (^), (^) and (|TT|), valid whenever the velocity 
is homogeneous, irrespective of the potential being or not being homogeneous, 
is 

E{t) = {v') = {(d^^f) = 2dt{^{x, t)). (21) 



^In this paper the symbol ~ is used for the leading-order term in asymptotic expansions, 
~ is used to relate two quantities which may differ by an order unity dimensionless constant 
and oc when the quantities are just proportional. 
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It relates the energy to the time-derivative of the mean potential, which does 
not dependent on the position x. 

We now recall some results about the Burgers equations which are rela- 
tively standard and will be used in what follows. 

When n > 0, the correlation invariant 

/CXD 
(f(x,t)t;(0,t))c/x = E(0,t) (22) 
-oo 

is time independent. Indeed, using the conservation form of the Burgers 
equation (^, it is checked that the time derivative of the integral in (|2^) , 
limited from —X to +X is expressed in terms of correlations between quan- 
tities evaluated at points a distance X apart. Such correlations vanish as 
X oo because, with our assumptions, the initial conditions have decreas- 
ing correlations at large separation and this property will be preserved.! Note 
that J is just the value of the energy spectrum at = 0. 

A straightforward consequence of the Burgers equation with a homoge- 
neous velocity field is 

dtB,{x,t) = ^d:,S^{x,t) + 2udlB,{x,t), (23) 

which relates the correlation function B^{x,t) and the third order structure 
function S^i^x) = {[v{x,t) — v{x,0)]^). This equation is the analogue of the 
Karman-Howarth equation for Navier-Stokes turbulence. Just as the latter 
implies Kolmogorov's four- fifths law (see, e.g., Frisch 1995), eq. ( p3D implies 
a "twelve law" : in the limit of vanishing viscosity and at scales much smaller 
than the integral scale, we have 

Ssix) = -12ex, (24) 

where e = dt{v'^/2) is the mean energy dissipation. 

Eq. (^) or (plD , in the limit of zero viscosity, puts a strong constraint 
on self-similar decaying solutions both for Navier-Stokes and Burgers tur- 
bulence. By self-similar decaying solutions, we understand the following : 
when the velocity v and the space coordinate x are rescaled by suitable 

■^This is expected because influence cannot propagate faster than the initial velocities. 
So, correlations over a distance X should not exceed a value oc e"^^/^^'^ /(to-i.) ) ^ ^ should 
be possible to prove this rigorously for the kind of initial conditions assumed here. 
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time-dependent scaling factors u{t) and L{t), the (single-time) statistical 
properties of the solution become independent of the time. This definition 
implies in particular the following scalings for the correlation function and 
the third-order structure function of the velocity : 

B^{x,t) = u\t)B^(^-^^ (25) 
Ss{x,t) = u\t)sJ-^], (26) 



with dimensionless time-independent functions and S3. Substituting (pSl)- 
( ^61) into ( p3|) for z/ — > and demanding that B^ and S'3 be time independent, 
we obtain two first-order differential equations relating u{t) and L{t), namely, 

u L ■ , . 

- = ai — , L = a2U, (27) 
u L 

where the dot denotes time differentiation and ai and 0:2 are two dimen- 
sionless undetermined constants. Note that exactly the same equations are 
obtained with a self-similarity ansatz for Navier-Stokes turbulence. Eq. (^) 
has power-law solutions^ but, as already pointed out by Kolmogorov (1941), 
the exponents in the power laws cannot be determined solely from the self- 
similarity ansatz. Additional ingredients are needed. 

For what follows it will be convenient to give a precise definition of the 
scaling factors. We choose u{t) to be the r.m.s. velocity 

u{t)^{v\xM'\ (28) 

and we define the "integral scale" L{t) by taking 02 = 1 in (pTj), that is 

L{t) = u{t). (29) 

Let us finally note that in Fourier space the self-similarity ansatz (|25|), 
together with (p9|), translates into 

E{k,t) = ^E{kL{t)), (30) 



''^In what follows, we shall be concerned with decaying solutions which satisfy the self- 
similarity conditions (p5|)-(p6|) and only in an asymptotic sense, for long times and to 
leading order. We shall then see that u(t) and L{t) may take the form of power laws with 
logarithmic corrections. 
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E{k) = — / B^(x) exp(ikx) dx, (31) 
27r J-oo 



where our choice of normahzation of the energy as u'^{t) imposes that / E{k) dk 
1. 



3 Phenomenology of the decay 

3.1 Permanence of large eddies and self-similar decay 

An important role in the decay at long times is played by the principle of 
permanence of large eddies (PLE). We shall give two formulations almost - 
but not quite - equivalent, one in physical space and the other one in Fourier 
space. 

PLE : physical space formulation. We assume that the Gaussian ini- 
tial velocity has a correlation function with a power-law behavior at large 
distances : 

~ a^C„|a;r"-\ |x| ^ oo, (32) 

with an exponent n > —1. Then, (^) remains true at any later time, with 
the same exponent n and the same constant a. 

PLE : Fourier space formulation. We assume that the Gaussian initial 
velocity has a spectrum with a power-law behavior at small wavenumbers : 

E(fc)~a2|A;r, |fc|-^0, (33) 



with an exponent 2 > n > — 1. Then, ([33| ) remains true at any later time, 
with the same exponent and constant. 

In its physical-space formulation PLE can be supported by the following 
argument. From (p^, it follows that the characteristic velocity associated 
with a large separation x is 

Uf){x) ~ a\x\ (34) 

so that the corresponding "turnover time" , the characteristic time for signif- 
icant changes through nonlinearity, is 

Uq\X} 
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which grows without bound at large separations. Hence, at any finite time, 
eddies of sufficiently large size should remain essentially as they were initially. 
This simple argument could be upset if two ffuid particles, initially strongly 
correlated (say, separated by not more than one initial integral scale), were to 
propagate this correlation to a large separation x in a finite time.i However, 
with Gaussian initial conditions this occurrence has a very small probability, 
which can be bounded by oc exp[— (l/2)a;^/(cr^t)^], as explained in the foot- 
note following the derivation of the constancy of J (eq. p^)) in the previous 
section. Thus, the initial correlation function at large separations should 
remain unaffected to leading order. 

In its Fourier space formulation, PLE has already been proven in the pre- 
vious section, for the case of the spectral index n = 0. When n < 2, PLE is 
a consequence of the observation that the time-rate of change of the energy 
spectrum, the so-called transfer function, is typically proportional to k"^ at 
small wavenumbers, because of the beating interactions of pairs of Fourier 
components with nearly opposite wavevectors. We may also derive the con- 
stancy of the coefficient of just from the constancy of the coefficient of 
at large x. 

There are reasons why the two formulations are not quite equivalent. 
First, the n = result in Fourier space is not an instance of PLE in physical 
space. Indeed, when n is an even integer, (^) implies a non-analytic spec- 
trum at small k with logarithmic corrections. Second, when n > 2 and the 
initial correlation function has the behavior (|3^) at large separations, the en- 
ergy spectrum at any later time will begin with an analytic (time-dependent) 
fc^ term with a subleading non-analytic correction proportional to \k\^ (with 
a time-independent coefficient). 

We now show how PLE, used as a "boundary condition" at large x or 
small k in connection with the self-similarity ansatz (^) and (pQ]), allows us 
to derive the law of temporal variation of the integral scale and the energy. 
We first give the argument in physical space. Let us assume that the distance 
beyond which PLE applies is within the scope of the self-similarity formula 
(P^l). It then immediately follows from (pSj), ( pSj ) and (^) that B{x) oc 

^For the case of incompressible Navier-Stokes turbulence, pressure effects can instan- 
taneously affect the large-scale correlations and thereby invalidate PLE, at least in the 
anisotropic case (Batchelor & Proudman 1956). 
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\x\ " and 

L{t) ~ (at)^ (36) 

and thus, by (pQ]), 

„ 4 2(n + l) 

E{t) = u^{t) a^t~^+^ . (37) 

In Fourier space the self-similarity ansatz (|30D, together with (|33| ) gives 
i?(fc) oc Ifcl'^, at small and the same relations for integral scale and the 
energy as written above. Clearly, this argument cannot be applied with 
initial data such that the spectral index n >2, since the later spectrum has 
now a k'^ dependence at small k with a time- dependent coefficient. 

Actually, we shall now show that the laws (0) and (|37|) apply only when 
n < 1. 

We begin with the case — 1 < n < 1 when the initial potential has homo- 
geneous increments. Many aspects of this case are well understood, thanks 
in particular to Burgers' own work (1974; see also Gurbatov et al. 1991). 
The phenomenology is quite simple. Increments Aipo[L) of the initial po- 
tential over a distance L = x — a can be estimated from the square root 
of the structure function Sq^{L) of the potential (|^). For a given position 
X, the maximum in (P)-(0) will come from those as such that the change in 
potential is comparable to the change in the parabolic term 

/sWW~y. (38) 

At large times L is also large and we can use (^) to obtain precisely (^) 
and, hence, the energy law (^7|). Alternatively, we could argue that when t is 
so large that the parabolas appearing in (|^ have a radius of curvature much 
larger than the typical radius of curvature of features in the initial potential, 
we can plausibly replace that initial potential by fractional Brownian motion 
of exponent h = (1 — n)/2 (in order to be consistent with (pO])). Without loss 
of generality we may assume that this fractional Brownian motion starts at 
the origin for x = 0. This function is then statistically invariant under the 
transformation x \x and ipo — > X^ipo- It is then elementary, using (P)-(0), 
to prove that a rescaling of the time is (statistically) equivalent to a suitable 
rescaling of x distances and of ip{x,t). This implies (|36|) and (^Tf).! 

^For details, see Vergassola et al. (1994, Section 4) and Avellaneda E (1995). 



13 



3.2 Elementary derivation of the Kida law of decay 

The case n > 1 when the initial potential is homogeneous is more difficult. 
Because of this homogeneity, the structure function becomes constant at large 
distances. If we then formally use we obtain a law for the energy 
and a t^^^ law for the integral scale. This is almost, but not quite, the right 
result. There are logarithmic corrections, ffist discovered by Kida (1979), 
which can actually be captured by phenomenology, as now explained. 

We shall estimate the energy in terms of the mean potential at x = 0, 
using ( pID and the maximum representation (|D-(^. We assume that, at 
long times t, the distribution of values of the potential is peaked around 
its (non-vanishing) mean value (ip) = {ip{0,t)). This quantity is thus the 
typical maximum value of ipo{a) —a?/{2t). Equivalently, we ask: what is the 
value {ip) by which we should vertically shift the parabola a? / {2t) so that the 
graph of ipoia) has a substantial probability of being everywhere below the 
shifted parabola ? To estimate this probability we divide the a-axis into equal 
intervals, /i, I2, I3, ■ ■ ■ , with a common length of about one initial integral 
scale Lq. The initial potential takes approximately independent values in the 
different intervals. Thus, the probability of no intersection with the shifted 
parabola is (approximately) the product of the probabilities qi, q2, . . .of no 
intersections in the various intervals : 

^no inters. ~ 91^2 • • • • (39) 

Each of the g^s is of the form 1 — pi, where pi is the probability of at least 
one intersection in the interval Jj. This can be estimated as 

Pi^W (^{^) + . (40) 

Here, W{iIj) is the single-point cumulative probability of ipo, integrated from 
ip to infinity and Oj is any point in the interval Jj, since the parabolic function 
changes very slowly with a when t is large. All the PiS are small and thus 

^no inters. ~ 1 - Pi - P2 • • • • (41) 

The sum of the PiS can be approximated by an integral from —oo to +00. 
When (ip) is taken very large this integral is very small. As we decrease 
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{ip) , the probability of having an intersection somewhere becomes significant 
when {ip) is chosen such that 

Since we assumed Gaussian initial conditions with r.m.s. value (T^,the 
cumulative probability W{ip) is approximately exp [— (l/2)('?/'/cr^)^]. The 
leading order of the integral in (^21) is obtained by keeping only the first two 
terms in the expansion of the square of {ip) + a? /{2t). Thereby, we obtain 

W~a,ln^/^(^), (43) 

where 

tnl = ^ (44) 

is the characteristic nonlinear time. Using (pT]), we then obtain the Kida 
(1979) log-corrected 1/t law for the energy decay: 

E(t)~rV^ln-^/2(^). (45) 

Finally, assuming (|^), we obtain Kida's (1979) result for the integral scale: 

m t^l\a,flH^-y^ {Ly (46) 

We must, of course, stress that this derivation is only heuristic. We shall 
come back to the issue of a more systematic derivation (Sections ^ and |^), 
including estimates of the subleading corrections. 

It is straightforward to extend this phenomenology to initial potentials 
which are homogeneous and non-Gaussian, as the only ingredient used is 
the p.d.f. of ipQ. For example, when the p.d.f. has a tail oc exp(— l?/;!^), the 
logarithmic correction to the 1/t law has the exponent (1 — 13)/ (3. Hence, 
when n > 1 the law of decay could be very sensitive to the functional form 
of the p.d.f. of the initial potential. Such heuristic considerations for non- 
Gaussian potentials have not yet been tested by more systematic theory and 
simulations. 
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Figure 1: Sketch of the energy spectrum at long times when 1 < n < 2. 



3.3 The Kida law and the issue of self-similarity when 

1 < n < 2 

One of the most questionable aspects of our phenomenological derivation is 
its validity when n is not an even integer. Indeed, this gives rise to long- 
range algebraically decreasing correlations, which are neglected in this phe- 
nomenology and which Kida (1979) explicitly ruled out. Actually, we shall 
see in Section |^ that the Kida law remains valid nevertheless as long as the 
initial potential is homogeneous, that is for any n > 1. 

We now address for the first time one of the central issues of this paper. 
The predictions (|36D and (^) for E{t) and L{t), obtained by assuming self- 
similarity and using PLE as a "boundary condition" at large x (or small k) are 



not consistent with the Kida law (45)-(46). Because of the time- dependence 
of the coefficient of k"^ in the spectrum at small k, the applicability of the 
former is obviously restricted to n < 2. Since the latter applies when n > 1, 
there may be a contradiction in the region of overlap 1 < n < 2. Actually, 
we shall see that n > 1 leads to a self-similar decay of the solution in an 
"inner" region of not too large distances (or, in Fourier space, of not too small 
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wavenumbers), controlled by the integral scale L{t) p6|). In this inner region 
correlations of ip decrease like a Gaussian function and the Fourier signature 
is a range of wavenumbers at which the energy spectrum is oc /c^. In addition, 
if n > 1 and is not an even integer, there exists an "outer" frozen region at 
which PLE applies. Thus, there is, at very large distances, an algebraic tail 
in the correlation function and the spectrum has a time-independent |/c|" 
contribution at small k. When 1 < n < 2 this is the leading term at small 
k so that the complete spectrum has three ranges: proportional to l/cl" at 
very small k, then proportional to A;^ at intermediate wavenumbers (up to 
approximately the inverse of the integral scale L{t)) and, finally, proportional 
to k~^ at high k (the usual signature of shocks), as shown in Fig. H We shall 
return to the issue of the global picture in Section |^. 

4 Homogeneous initial potential : previous 
work 

In this section we discuss previous work dealing with the limit of vanishing 
viscosity for initial potentials which are homogeneous and Gaussian. The 
main papers under discussion are Kida (1979), Gurbatov & Saichev (1981; 
see also Gurbatov et al. 1991), Fournier & Frisch (1983) and Molchanov, 
Surgailis & Woyczynski (1995). For initial conditions such that the velocity 
field is homogeneous but the potential is not, or such that the initial condi- 
tions are not Gaussian, the reader is referred to Burgers (1974), Albeverio, 
Molchanov & Surgailis et al. (1994), and Esipov & Newman (1993), Esipov 
(1994), Hu & Woyczynski (1996) and references therein. 

Let us first state the essence of what is at presently known to hold with 
certainty. In the limit of vanishing viscosity, as the time t tends to infinity, 
the statistical solution becomes self-similar in the sense of Section 0. The 
integral scale L{t) and the energy E{t) are given, to leading order, by 

L{t) ^ (47) 

Fit) ^''^^^^''^'{^)^ (48) 
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where 

tn. - ^ = ^, L, ^ ^. (49) 

The non-dimensionahzed self-similar correlation function B^[x), defined by 
(|25D , which is a function of x = a;/L(t), is given by 

B,{x) = ^ (xP(x)) , (50) 



where 



1 dz 
^""^ = 2 loo ^ (^) exp [(^±^] + g (^) ^^^^ 



exp 



and 



„2 



^(^) = / e-~ rfs. (52) 



oo 



Note that the properties of the self-similar state are universal in so far as 
they are expressed solely in terms of two integral characteristics of the initial 
spectrum, namely the initial r.m.s. potential and r.m.s. velocity a^. I It 
may be shown that the function P{x) is the probability of having no shock 
within an Eulerian interval of length xL{t) (Gurbatov, Malakhov & Saichev 
1991). 

Important consequences of (|50D-(pTD are the asymptotic behavior of the 
non-dimensionalized correlation function By[x) at large x and of its Fourier 
transform, the non-dimensionalized spectrum E{k) ( PT| ) at small k = kL{t) 
(Gurbatov, Malakhov & Saichev 1991) : 



Bv{x) ~ xexp 1^-- , |x| ^ oo, (53) 

Eik) ~ /i^P, 1^1 ^0, (54) 
1 

= _ / x^p(x)dx ^ 1.08. (55) 
vr Jo 

Most of these results have already been obtained by Kida (1979). He was 
interested in the long-time behavior when the initial potential is homogeneous 



^Observe that the spectral exponent does not directly enter, in contrast to what hap- 
pens when n < 1 (cf. (|3^)-(37)). 
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and the initial correlations decrease rapidly at large separations, without 
particularly emphasizing the Gaussian assumption. He introduced a model 
of discrete independent maxima whose p.d.f. had three free parameters (an 
exponent and two constants); their relation to the properties of the initial 
conditions (say, the spectrum) were left unspecified. For the case of a p.d.f. 
with a Gaussian tail, he obtained the functional form of the results given 
above. 

At long times, the contacts of the parabolas, in the maximum repre- 
sentation (§)-(0) of the solution, are near large- value local maxima of the 
initial potential. Kida's assumption of independence of the successive max- 
ima implicitly assumes a Poisson process. Gurbatov & Saichev (1981) who 
conjectured the asymptotic existence of a Poisson process! and Molchanov, 
Surgailis & Woyczynski (1995) who proved it, showed that, in the x-%Ij plane, 
the density of the points is uniform in the x-direction and exponential in the 
^/'-direction. This permits the calculation of the one- and two-point p.d.f. 's 
of the velocity (Gurbatov & Saichev 1981). Molchanov, Surgailis & Woy- 
czynski (1995) also calculated the full iV-point multiple time distributions. 
Their proof made crucial use of some fine properties of extremes of Gaussian 
processes (Leadbetter, Lindgren & Rootzen 1983). i From the point of view 
of the present paper it is important to stress that Molchanov, Surgailis & 
Woyczynski's (1995) result is valid for any n> 1 and not just for n > 2.li2l 

Assuming a Poisson distribution, Gurbatov & Saichev (1981) and Fournier 
& Frisch (1983) showed that the statistical properties of the points of con- 
tact between the parabolas and the initial potential can be obtained from 
the statistical properties of their intersections, whose mean number can be 
calculated using the formula of Rice (1954; see also Leadbetter, Lindgren & 
Rootzen 1983; Papoulis 1991). Thus, they could express the parameters in 
the asymptotic formulas in terms of the r.m.s. initial potential and velocity. 

None of the above theories give a control over how fast in time the asymp- 
totic results are reached. Furthermore, the rescaling involved in these theories 
can only capture what we have called in the previous section the "inner re- 



®See also Fournier & Frisch (1983). 

^What they proved is not exactly in the form of the long-time asymptotic relations given 
above but an equivalent result in which the time and space variables and the amplitude 
of the initial conditions are rescaled. 

^°They assume that the initial correlations of V' decrease faster than the inverse of In |x| 
at large |a;| and can be Taylor-expanded to fourth order at small |a;|. 
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gion", where PLE is irrelevant. We shall show in the next section that, for 
finite long times, the Poisson distribution is just an approximation which can 
be systematically refined, so as to address such issues. 



5 Systematic derivation of the Kida law for 
the energy decay using the Bahan— Schaeffer 
formula 

The (mean) energy is given by (|2l|) in terms of the time derivative of the mean 
potential (ip). Here, we use the method of Fournier & Frisch (1983) and the 
Balian-Schaeffer formula (|A.18|) , established in Appendix ^ to derive the 



asymptotic behavior of the mean potential for the case of a homogeneous 
initial potential, that is, n > 1. 

Let P{A,t) denote the probability density of at time t and let Q{A,t) 
denote the (cumulative) probability to have ip < A, given by 

A 

Q{A,t) = J P{A',t)dA'. (56) 

— oo 

The mean value of the potential at time t is then expressible as 

oo 

{^) = J AP{A,t)dA. (57) 

— oo 

An integration by parts gives 

oo 

{^) = /[I - Q{A t)] dA- J Q{A, t) dA. (58) 



We shall see that the second integral in ( |58D may be bounded by a time- 
independent constant, while the first grows; so it is only the first integral 
that will be of importance in what follows. 

From (|^)-(0) and the homogeneity of the initial potential, we have 



Q{A,t) = Prob 



My) - 1^ - ^ < 0' all y 



(59) 
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The argument of the probabihty in (|59| ) can be written as the condition 
of non-intersection of the following two graphs, shown in Fig. |]: 

Gi: ^poiy), 

G2: y^A + ^^. (60) 

We can immediately use this geometric picture to eliminate the second 
integral in (^) from further discussion. Indeed, the probability that Gi and 
G2 never intersect is obviously less than the probability that ipo{y) is less 
than A at the single point y equal to zero. The integral over A from minus 
infinity to zero of this quantity is finite and bounded by a quantity of the 
order of a^. This follows from the observation that, for negative A, the 
probability Q{A,t) given by (^) is less than the same expression evaluated 
at the single point y = 0. Hence, the contribution of the second integral in 
( |55D is bounded by a time-independent constant and will not contribute to 
the logarithmically growing leading-order mean potential. 

We now treat the intersections of the graphs Gi and G2 as "particles" 
and apply the Balian-Schaeffer formula ( |A.18| ) of Appendix ^ to calculate 
the probability of having zero particles (no intersections at all). We shall 
limit ourselves to the case where Gi intersects G2 going upward (see the 
footnote on p. 30 for what happens when all intersections are counted.). 
Such intersections are depicted by circles in Fig. |^. 

From (H), (13) and (^1|), we have 

/ 00 t \n \ 

Q{A,t)=po = exp[J2^-^cJ, (61) 

\n=l / 

where the cumulants c„ are constructed from the n-particle distribution func- 
tions Fn as explained in Appendix 0. 

In particular, Ci = (A^) = / Fi{y)dy, the mean number of particles on the 
whole real line. We now claim that we have the Rice (1954) formula: 

{N) = J F,{y)dy=(^J 5 ($ (y, t)) , (62) 

where $(?/, t) = ipoly) — ^ — A and the shorthand denotes the positive 
part of X, i.e. equal to a; if x is greater than zero and zero otherwise and 5 is 
the delta-function. 
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Indeed, the integral in (|62|) only gets a contribution when the function $ 
vanishes. This happens precisely when the two graphs Gi and G2 intersect. 
Let y* be a point of intersection. An integral over a small interval around 
the intersection then contributes to (|6^) a quantity \^y{y*,t)\~^{^y{y*,t)}~^ 
. If the derivative of $ at ?/ is positive, that is, if Gi intersects G2 going 
upward, then the contribution is one. If, on the other hand, Gi intersects G2 
going downward, then the contribution is zero, because {^'y{y*, t)}^ vanishes. 
The integral in (|62D therefore counts the number of upward intersections 
(upcrossings) between the two graphs. Note that the assumed smoothness 
of the process ipo{y) (which follows from the fast decrease at infinity of its 
spectrum) renders the applicability of the Rice formula robust. 

The integral / Fi{y)dy is thus determined by the joint probability density 
W^^^i (^$, y^ of $ and its derivative $' at the point y : 

{N) = J F,{y)dy = Jdy J d<^' {<^'}+W^ ^> ($, y) . (63) 

We have similar expressions for the functions F„ (?/i, . . . , ?/„). Let us in- 
troduce the n-point joint probability density W^^^ {ipiy • • • > i^n] iJi, ■ ■ ■ , Vn) of 
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the initial potential ip and the initial velocity v, the latter being equal to 
minus the space derivative of the initial potential. 
Then we can write, analogously to (^) : 



Fn {yi, ...,yn)=J dvi... dVnW^^y ^ 

■{vi - yi/ty{v2 - y2lty ■■■{vn- yn/t} 



^ + A,...,y^ + A;vi,...,Vn 



(64) 

Since the initial potential is Gaussian, all the joint probability densities 
are completely determined by the initial two-point correlation function : 



Bq^ (x) = (^0 {y + x) ipo (y)) 



(65) 



In our case we need the joint probability densities for (initial) fields ip at 
various points and their (negative) space derivatives v at those same points. 
They are given by 



W^^v {ipl, . . . ,1pn;Vi, . . .,Vn) 



(27r)^" det M 



1/2 



• exp 



where M is the enlarged correlation matrix 



Bq^ {yt - ye 



BL {ye 



B'L (ye -ye') = Bo^ {y^ - yp,) 



In the last line we have introduced the initial velocity correlation Bq^. 
The diagonal elements of this matrix are 



a. 



(66) 



(67) 



(68) 



Here, and are the variances of the initial potential and the initial 
velocity, respectively, as previously defined. 

We now return to the Balian-Schaeffer expression (|6lD for the cumulative 
probability of the potential. Inspection of (|58|) shows that the mean potential 
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comes predominantly from those arguments A such that Q{A,t) is not too 
close to unity or, by (|6l|) , such that the alternating sum of cumulants is not 
much smaller than unity. As A grows, the changeover of the alternating sum 
to a value of order unity happens around a value which we shall denote A* (t) . 



5.1 The Poisson approximation 



We are interested in the energy decay at large times and it will be shown 
that in this case the first cumulant Ci = /i is the leading term in (|6T|) and 
that, to leading order, we have the Poisson approximation: 



Q{A,t) ^ exp{-ci} 



(69) 



It is so called for the natural reason that if the upcrossings obeyed a Poisson 
distribution, then (|69|) would be true exactly. The rest of this subsection 
contains a computation of the first cumulant ci. 

The function Fi is determined by the single-point probability density of 
the initial potential ip and the velocity v, the subscript being dropped for 
conciseness. For a homogeneous random function, these two variables are 
uncorrelated. We then obtain from ( |63|) and (16^) 



Fi{yi 



exp 



iyl/2t + Af , 



at, 



{v-^}+dv. (70) 



It is now useful to introduce the auxiliary function 



Z(x) 



1 \ ^/^ 

-j l{v-x}^e- 



2 dv. 



(71) 



which can obviously be expressed in terms of the error function. We only 
need to use that, for x — > 0, the function tends to Z{0) = (27r) ' and that, 
for X — * oo, we have Z{x) ^ x. 

In terms of Z{x) the mean number of upcrossings is given by 



Ci = y Fi (yi) dyi = J dyi Z (^j^^ y^l/27ro-| exp 



{y!/2t + Af 



(72) 

The cumulant ci is a dimensionless number. It is therefore preferable to 
rewrite ( |72D in a manifestly dimensionless form. By (0), o^ja^ is equal to 
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Lq, the initial integral scale. In (^) we introduce the new dimensionless 
variables zi = Ui/Lq and a = A/a^ and the dimensionless nonlinear time 
tni = Lo/a^ already defined (^^j). Eq. (^) can then be written 




Ci = / dziZ ( I -i/l/27rexp 



,2 n2 

i-^ + a 



2(tAnl) 



(73) 



It is seen that (|73|) involves two spatial scales : when zi is small compared 
to the "prefactor scale" t/t^i we can simply substitute Z{0) for Z{zi/{t/tn\)); 



when zi is large compared to the "Gaussian scale" l{a, t) ~ y (t/tni) / o there is 
a large damping by the term in the exponent quadratic in Zi. If we assume, 
for the moment, that the typical value of a* = A* /a^ where ci becomes 

order unity does not grow in time faster than '^t/t^x^ it follows that it is 
self-consistent to assume that the prefactor scale is much larger than the 
Gaussian scale, so that we may safely use Z{Q) in (|73D. 

The remaining integral in (^) cannot be carried out explicitly, but it 
is nevertheless easy to evaluate it to relevant leading order. We have two 
terms in the exponent that depend on Zi : one quadratic and one quartic. 
The quartic term is larger for Zi much beyond the Gaussian scale, but then 
the damping is strong anyway; so, the correction from the quartic term is 
small, of relative order and thus negligible. It is then easily shown 
that, for those large values of a ~ a* which give the dominant contribution 
to the mean potential at large-time, we can use the following leading-order 
approximation for the cumulant Ci : 



ci(At) ^ ( r^V^'e-"'/', A = a^a. (74) 



aV2 



We now determine a*{t). For this, we assume that Ci is the dominant term in 
the Balian-Schaeffer formula (|6T|) for the probability of no intersections, as 
will be shown in Section ^.2[ a*{t) is then determined by the condition that 
ci{a^a*,t) be order unity. To avoid undetermined constants, we define a*{t) 
by the condition ci{a^a*,t) = 1, that is through the transcendental equation 

+ \ 1/2 

* ^ e-('^*)'/2 = 1, (75) 



tni27ra* 



11 



Actually, it will be seen to grow only as a square root of the logarithm of the time. 
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from which it is easily shown that, to leading order, at large times 

We now return to the cumulative probability of Q{A,t). Using the Pois- 
son approximation ( |69D and the leading-order approximation (^) for ci, we 
obtain the following doubly exponential distribution :0 

Q(A,t)^exp{-exp{-C}}, (77) 
where ( is obtained from a by shifting and rescaling : 

a = a-(l + ^), (78) 

Since a*{t) is logarithmically large at large times, it follows from (|7^-(0) 
that the values of a = A/a^ contributing to the mean potential ( ^Sf ) are 
concentrated near a*{t) with a small relative width ^ _ Aa i 



2 



A* a* (a*) 

In 2^1^. Hence, to leading order, the mean potential is given in the Poisson 



_t 

TTt 

approximation by : 



(V;(t))^a^a*^a^W/2(^), ^^g^ 

which implies the energy decay law (|48|) . 

Notice that, within the Poisson approximation, the energy decay is fully 
determined by the local properties of the initial potential correlation function 
Bq^ (x), namely, cx^ = Bq^ (0) and = —Bq^ (0) and does not involve its 
large-scale behavior. 



5.2 Corrections to the Poisson approximation 

So far we have only reproduced a known result. The Balian-Schaeffer formula 
allows us to find the applicability condition for the Poisson approximation 
by estimating contributions of higher-order cumulants to (^TJ). 

^^This doubly exponential distribution is equivalent to the distribution of maxima be- 
ing Poisson with uniform density in space and exponentially decaying in a (Molchanov, 
Surgailis & Woyczynski 1995). 
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We thus turn to the calculation of the second cumulant C2 appearing in 
(|6lD , given in terms of the one- and two-particle densities by (|A.17|) , repeated 
hereafter for convenience : 



C2 = / / [^2(2/1, 2/2) - Fi(yi)Fi(y2)] dy^dy^. (80) 

The function F2 is now given by (Q). Our purpose in this section will only 
be to show that the contribution of C2 - and of the higher-order c„s - are 
subdominant compared to that of ci. For this purpose, it will be enough 
to estimate the order of magnitude of C2 (for values of a ~ a*{t)). Thus, 
there is no need to derive the full leading-order expression for C2 with all the 
numerical constants. 

From (|64D-(|67D it follows that, unlike Ci, the second cumulant C2 involves 
the entire initial correlation function Bq^{x). Let us denote by Acorr a typical 
correlation length for the initial potential.EJ We shall now distinguish those 
contributions to (|80D where the two arguments yi and 1/2 are separated by 
not more than Acorr and those where the separation is larger or, possibly, 
much larger, a case of interest when the initial correlation function has a slow 
algebraic decrease. The former will be called "local" (loc.) and the latter 
"non-local" (nloc). 

Using (|6^) and (^O]) we can estimate 

- j j F2 (z/i, 1/2) dy^dy2 ~ Acorr / {Fi{yi)f dy,. (81) 
Introducing as in (|73|) new (dimensionless) variables zi and a, we get 

The integral in (|8^) is practically the same as for the first cumulant, so we 
have 

c'riAt) ~ ^ (t^T^--'^ a = (83) 

The main differences with the expression ([7^ ) for ci is (i) the presence of a 
factor Acorr / Lq and (ii) the fact that the argument of the exponential is now 

""^^We could have taken Acorr = Lq but keeping the distinction has some interest, as we 
shall see. 
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—a^ instead of — From ([75| ) and ( ]83| ) we find that, around the value 
a*{t) which give the dominant contribution to the mean potential, the second 
order local cumulant is roughly 



Joe. / ^corr 



Lo J V^nia*(t), 



^1/2 

(84) 



It easily follows that the local second order cumulant may be neglected com- 
pared to Ci as soon as 



A 



2 



. (85) 

Here, we make a short digression. There are instances where (Acorr/-^o)^ 
can be large. Consider for instance a quasimonochromatic signal with a 
center wavenumber ko ~ Lq^ and a width Ak ~ [Acorr]^^ ^ ^o- In this 
case, the condition t/tni ^ 1 is not enough for the Poisson approximation 
( |5UD to hold. The asymptotic regime described by Kida's law is then obtained 
much later. As long as tni -C t -C tni(Acorr/-f'o)^) the energy decays much 
faster: E{t) ~ Lg/t^ (Gurbatov & Malakhov 1977). The physical reason for 
this is a strong correlation of the shocks in the early stage of the evolution, 
which prevents the rapid merging of shocks.lli 

Let us now turn to the non-local contribution to C2 which is important 
only when n is not an even integer, so that there is an algebraic tail in the 
correlation function : 

/A 

Bo.^{x)r^al^-^j , x>Acorr. (86) 

Owing to this, when 1 < n < 2 the main contribution to C2 comes from the 
non-local region where \yi — 1/2I 3> Lq while \yi\, I1/2I < L{t). Here, L(t) is 
the integral scale at time t, given by (^61). 

The main technical difficulty now is that we have to express the inverse 
matrix of correlations that appears in the exponent in the Gaussian joint 
probability densities (|6^)-(|67|). Actually, there is no need to fully invert the 
matrix. It is sufficient to note that when the points are far apart, correlations 
become weak and the off-diagonal elements will be small compared to the 



Similar remarks can be made when the initial energy spectrum has a very long plateau 
at intermediate wavenumbers. 
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diagonal ones. We may therefore invert perturbatively to first order in the 
off-diagonal elements. Furthermore, the velocity-velocity correlations and 
the velocity-potential correlations decay much faster with separation than 
the potential-potential correlations. Therefore, it is consistent to keep in the 
inverted matrix only off-diagonal elements proportional to -BoV'(z/£ ~ Vi')- -^^ 
before we can substitute for the function Z(x) its value at the origin, (27r)~^. 
After all these simplifications we can write the following expression for the 
non-local part of the second cumulant C2 : 



nloc. 




dyi dy2 ■ exp 

Bq4, {yi - y2) 



1 



a. 



(87) 

Here, the subscript NL (non-local) is a reminder that we only integrate over 
points yi and y2 sufficiently far apart. The next step is to expand the second 
exponential, which contains the correlation function, to linear order. As in 
the evaluation of Ci, we change the variables yi, 2/2 and A into the dimen- 
sionless variables zi, Z2 and a, respectively and obtain, using (R^) : 



nloc. 



A 



corr 



n-l 



[ dzi dz2 ■ exp < — f — r 
Jnl ^ I 2 V2(tAni) 



-|- a 



2{t/t 



nlJ 



-|- a 



Fi - Z2 



1-n 



2(t/t 



+ a 



nl 



2{t/t 



nl 



(88) 



the integration being now over the region such that — 2;2| > Acorr/-^o- 



The integral (|88| ) is estimated in the following way: the exponentials 
effectively delimit the region of integration to a circle in the (zi, Z2)-plane 
with radius l(t) ~ (t/(atni))^''^- Inside this circle, most of the factors in (^) 
are of order unity. Therefore, we can estimate the non-local part of C2 as 



nloc. 



A 



corr 



NL 



dzi dz2 ■ e 



d \Z\ 



Z2\ , 



(89) 
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where NL is now the region 



t 



\Z2\ < 



t 



A 



corr 



(90) 



When n > 2 the non-local component of C2 is much smaller than the local 



one and we may ignore it. When 1 < n < 2 the integral in (39) is readily 
evaluated, and we obtain the estimate 



nloc. 



A 



corr 



Finally, estimating c^^'-'^- for A ~ A*{t), we obtain 



2 — n 
2 



(91) 



nloc. 



A 



corr 



n-1 



t 



(92) 



Comparison with (plf ) indicates that the non-local part of C2 is much larger 
than the local one, but still small compared to unity, so that we may indeed 
neglect C2 compared to Ci.El Let us observe that when n is just slightly larger 
than one, an exceedingly long time may be needed before the asymptotic Kida 
law applies since C2 is then just marginally smaller than Ci. 

Using similar procedures we can show that cjj^*^' and c^*^^' for m > 2 
are even smaller than the C2-corrections. This establishes the validity of the 
Kida law (|^) at long times. 



6 The correlation function at long times and 
the preservation of long-range correlations 

In this section we calculate the (spatial) correlation function of the poten- 
tial at long times for the solution to the Burgers equation, using a Balian- 
Schaeffer formula, adapted now to the evaluation of two-point quantities. 

""^^This statement holds provided only upcrossings or only downcrossings between the 
graphs Gi and G2 are counted. Otherwise, two things happen : (i) the mean number of 
crossings and hence ci is doubled; (ii) strong correlations are present between successive up 
and down large-ordinate crossings; these invalidate the Poisson approximation; actually, 
the coefficient C2 becomes large and exactly cancels the aforementioned factor two. 
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As we shall see, in evaluating the far tail of the correlation function, the 
dominant contribution will come from the second rather than from the first 
cumulant. 

For reasons of symmetry we take the following definition for the correla- 
tion function of the potential : 

B^[x, t) = (V'(x/2, t)^(-x/2, t)) - (^(x, t))\ (93) 

from which the the velocity correlation function is obtained as 

(94) 

To evaluate (|93|) we need, in addition to the one point probability density 
Pi(y4,t), the joint probability density P2{A, B]x,t) of the potential at the 
same time t for two points with spatial separation x. From (|93|), we obtain 

B^ = j J AB{P2{A,B;x,t) - Pi{A,t)Pi{B,t))dAdB. (95) 

We also introduce the cumulative probability Q2{A, B; x, t), related to P2{A, B; : 
by 

q2 

P2iAB;x,t) = ^^Q,(A,B;x,t). (96) 
Using this in and integrating by parts twice, we have 

B4x,t)= [Q2{A,B-x,t)-Qi{A,t)Qi{B,t)]dAdB. (97) 

J J — oo 

By the maximum representation solution (|^) we can express Q2 in the 
form 

Q2{A, B] x,t) = Prob (Gi never intersects G2) , (98) 
where Gi and G2 are now the following graphs, shown in Fig. |^: 



Gi: ipo{y), 

'{y-lf , , (y + if 



G2: y^ g{y) = min 



+ A, + B 



2t ' 2t 



(99) 
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B+(y-x/2)72t 




Figure 3: Intersections of the initial potential ipo with two parabolas. 



Thus, the procedure of calculation of Q2 is similar to that used in the pre- 
ceding section for calculating the mean potential, but with a more complex 
function g{y). Hence, we may again use the Balian-Schaeffer formula and 
obtain 



Q2{A, B- X, t) = exp X] ^-c„ 



where the c„s depend now on A, 5, x and t. 

The first cumulant ci{A^B) in the expansion ( |10CI| ) is now 



(100) 



cM.B)= Z{g'iy)/a.,) 



(27r(T^) 



1/2 



exp 



1 aivi) \ , 
— r "2/1- 



2 2a. 



(101) 



From (0) and Fig. ^ we see that the integral in ( |101|) can be divided into 
two integrals over pieces of the two parabolas intersecting at y* : 



ci{A,B) 



27lLn 



exp < 



(y+i 



2t 



B 



2 ^ 



> dy 
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+ 



2i 



+ A 



exp < 



> dy 



A-B 



y* 



2x 



-t. 



(102) 



(103) 



Here, we have assumed from the start that t ^ tni, to replace Z{x) by 
Z{0) = (27r)-V2. 

It may now be checked that ( |102| ) behaves in different ways depending on 
how the widths of the Gaussians compare to the separation x. 

As long as x is smaller than or comparable to the integral scale L{t) given 
by (13), the leading-order term in (|100|) is found to be Ci, that is, the Poisson 
approximation remains valid. One then derives the expressions (|50D-(p]) for 
the correlation function previously obtained by Kida (1979), Gurbatov & 
Saichev (1981) and Molchanov, Surgailis & Woyczynski (1995). Actually, to 
derive the self-similarly evolving correlation function given by (pO|), the only 
possible spatial scaling factor is (up to a multiplicative constant) the integral 
scale L{t) given by {^7\j. It may also be checked that L{t) is (up to a constant 
factor) equal to {tP^t))'/^ {v\t)y/^ 

When |x| ^ L(t) the situation is quite different. Now, the two parabolas 
intersect at a very large ordinate so that, to leading order, we have two 
independent integrations over two different almost complete parabolas. Thus, 
to leading order, we have 



Ci{A,B;x,t) ~ ci{A,t) + ci{B,t), x > L{t). 



(104) 



This, however, implies 



Q2{A,B;x,t) - Qi{A,t)Qi{B,t)c^exp{-ci{A,t)-ci{B,t)} 

- exp{-ci(A,t)}exp{-ci(fi,t)} = 0, (105) 

which, in turn, implies the vanishing of B^{x,t). It follows that the leading- 
order contribution cannot come from the Poisson approximation. 

This leading-order contribution actually comes from the second order 
cumulant 

C2 = J J[F2{yuy2) - F,{y,)F,{y2)]dy,dy2, (106) 
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Figure 4: The different regions of integration appearing in ([1061) . Near the 
diagonal yi = 1/2 we have regions I and I' of ( |107| ). Away from the diagonal 
we have regions II and 11'. 



where Fi and F2 are now the one- and two-particle densities constructed from 
the crossings of the Gi and G2 graphs defined in (|99D . It may be checked 
that when 3> L{t) the significant regions in ( |106| ) are: 



II : 



X 



y2 



ivi ~ 2' 



--) 



, X X ^ 

I : 2' ^2 ~ ~2' 



11': (1/1 



X 



X, 



(107) 



"2' ^ 



also illustrated in Fig. ^. 

In the regions I and I' we have the same estimates (|8^ -(|92D for the value of 
C2 as in the calculation of the mean potential. The corresponding contribution 
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to the correlation function goes to zero with separation even faster than the 
Ci-contribution and is not relevant. 

In the regions II and 11', for t ^ t^i and |x| ^ L(t) we obtain after 
integration over vi, V2, the following expression 



= — 2^ 



exp < JJq^ [X ) 



v 



4 



Ci(A)c,(5)(exp|^i^|j^|-l^^ (108) 



Just as in Section the joint probability density of A, B is localized in a 
narrow region near A ~ S ~ (^) ~ a^(lnt/27rtni)i/2_ ^Ve introduce in ( 110^ ) 
the new variables C, and ^ through 

A = (^) + |yC, ^ = (^) + , (109) 

and we observe that, when C, and are 0(1), we have, by (^) and (|79D , 

Ci{A) ^ e-^, ci{B) ^ e-«. (110) 

We now evaluate the correlation of the potential ( P7| ) for |a;| ^ -^(^)- For 
this we use ( |100| ), truncated to its first two terms ( |108| ) and ( |110| ), and the 
fact that Bq^{x) is small, so that we can expand exponentials with Bq^{x) 
in their arguments. We finally obtain : 

S^(x,t) ~ 5o^(x), t>tnl, |x|>L(t). (Ill) 

This establishes the preservation of long-range correlations, an instance 
of the permanence of large eddies (PLE), as defined in Section!^. 



7 The global picture at long times 

We are now in a position to derive the main results of this paper, some of 
which have already been presented in a phenomenological way in Section ^. 
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Our derivation of the Kida law ( P8| ) for the energy decay in Section ^ 
using the Bahan-Schaeffer formula, proves the validity of this law for any 
n > 1 and not just n > 2 as thought previously. Similarly, for any n > 1 the 
integral scale L(t) is given by p7|). 

More novel results emerge when considering the global structure, at long 
times, for the correlation function (in physical space) and for its Fourier 
transform, the energy spectrum. We shall consider only the situation at 
scales larger than the integral scale L{t). Indeed, we have already stated in 
Section ^that, at scales small or comparable to the integral scale, we recover 
known results (see Kida 1979; Gurbatov, Malakhov & Saichev 1991). 

We observe that for any n > 1 which is not an even integer the correlation 
function By{x) has two different asymptotic regions. In an "outer asymptotic 
region" at extremely large x, permanence of large eddies (in the form ( |1 1 Ij ) 
established in Section ^ implies that 

B^{x, t) ~ 5^(x, 0) ~ |x| oo. (112) 

There is also an "inner asymptotic region" where the correlation function has 
the (approximately) Gaussian decrease given by (|53D in dimensionless form 
with X = x/L{t). Using (pSf ) we can rewrite this in dimensioned form as 

which can be reexpressed in terms solely of L(t), using u(t) = L(t). The 
separation ls{t) at which the correlation function switches from its inner to 
its outer form, called here the switching distance, is obtained by equating 
( [ml) and dmD, yieldinjli 



/.(t)~L(t)((n-l)ln(^))'^'. (114) 

Hence, the switching distance is only logarithmically larger than the integral 
scale. c-3 Still, the fact that L{t) and ls{t) have (slightly) different scaling laws 



""^^To truly characterize the transition from the inner to outer asymptotic expansion it is 
necessary to perform matched asymptotic expansions which require more than the leading 
order. Equating the leading-order inner and outer terms correctly predicts the order of 
magnitude of the scale at which this transition takes place. 

^^When n is an even integer the algebraic outer region disappears altogether. 
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in time implies that, globally, the correlation function does not evolve in a 
self-similar fashion. 

This phenomenon becomes much more prominent if we consider the en- 
ergy spectrum, the Fourier transform of the correlation function. In dimen- 
sionless inner-region variables, the energy spectrum has the following form 
(Gurbatov, Malakhov & Saichev 1991) : 



E{k) ~ <^ ^ ^ , k = kL{t). (115) 




The k~'^ region is the signature of shocks and will not interest us further, 
while the k"^ region comes about because the (inner) correlation function of 
the potential has a non-vanishing integral from — oo to +00. In dimensioned 
variables the small-A; inner-region behavior of the spectrum is thus 

E{k, t) ~ k^ ^ A{t) k"^, kL{t) < 1, inner region, (116) 

t^ 



A{t) = alLl{^y\n-'/'(-^]. (117) 



where 

A(t^ = cT^T''J , ^ 

^tn\J \27rtni. 

Instead of a Gaussian region (in physical space), we have a spectrum with 
an algebraic k'^ region and a time-dependent coefficient A{t). 

We must now distinguish two cases. When n > 2, the k"^ contribution 
( |116| ) dominates everywhere over the \k\"' contribution, the latter being then 
only a subdominant term in the small-fc expansion.@ The more interesting 
case is 1 < n < 2, already sketched in Fig. |T[ The permanence of large eddies 
implies now that, at extremely small fc, 

£■(/;;, t) ~ a^l/cp, for A; 0, outer region. (118) 

This relation holds only in an outer region \k\ <^ ks{t) where ( |118| ) dominates 
over ( |116| ). The switching wavenumber ksit), obtained by equating ( |116| ) and 
( |118|) , is given by 



^^This subdominant term is easily seen in numerical simulations at very short times 
and moderately small wavenumbers, for example when n = 3 (A. NouUez, private 
communication) . 
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(120) 



Some comments are now made. First, it may seem paradoxical that the 
switching wavenumber ks{t) is not the inverse of the switching distance ls{t). 
The reason is that in physical space we have a Gaussian competing with 
a power law, whereas in Fourier space we have two power laws competing; 
this renders inapplicable naive Fourier phenomenology in which distances are 
roughly inverse wavenumbers. Second, let us define an energy wavenumber 
kLit) = L-^it), which is roughly the wavenumber around which most of 
the kinetic energy resides. From (^Tj) /^^(t) ~ (ta^)'^^"^ (ignoring logarithmic 
corrections). We then have from ( |119| ), still ignoring logarithmic corrections : 

Ks[tj / t \^2{2~n) 

kjt) ^ 

Hence, the switching wavenumber goes to zero much faster than the energy 
wavenumber : in Fourier space the separation of the two regions is much more 
manifest. This is why the results of the numerical simulations in Section || 
will be presented in Fourier space. 

Let us also observe that the ratio of the energy in the outer region to 
the total energy, a measure of how well the Kida law is satisfied, is equal to 

3(n-l) 

(t/tni) ""^ (up to logarithms) and thus becomes very small when t ^ tni, 
unless n is very close to unity. 

Thus, we have established the central claim of this paper, namely that for 
1 < n < 2 there is no globally self-similar evolution of the energy spectrum. 
Of course, as — > 2 the inner /c^ region overwhelms the outer region and 
as n — > 1 the converse happens, so that in both instances global self-similarity 
tends to be reestablished. 



8 Numerical experiments 

The purpose of this section is to obtain numerical confirmation of the rather 
unusual long-time regime predicted by the theory when the spectral exponent 
n is between one and two. There should then be three different power-law 
regions in the energy spectrum : 

E{k)oc\k\-^ for > L-^(t); 

E{k) oc \k\^ for ks{t) < \k\ < L-\t)] (121) 
E{k) oc for \k\ < k,{t), 
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where ks{t) and L{t) are given by ( |119| ) and (pT]), respectively. 

Clearly, such simulations require extremely high spatial resolution since 
three power-law ranges are present. Fortunately, for the Burgers equation 
there is no need to numerically integrate (|l|) : advantage can be taken of the 
maximum representation to directly construct, from the initial data, the 
solution in the limit of vanishing viscosity at any given time t. The principles 
of the method have been presented in detail elsewhere (She, Aurell & Frisch 
1992; Noullez & Vergassola 1994; Vergassola et al. 1994; Aurell, Gurbatov 
& Simdyankin 1996) and will only be briefly recalled here. 

The method exploits two observations. First, the maximum representa- 
tion (^) becomes a Legendre transformation after expansion of the quadratic 
term. Second, when the problem is discretized on a grid of N points, the 
search for the maximizing Lagrangian point a(x, t) can be done by a di- 
chotomic procedure which uses the non- decreasing property of the inverse 
Lagrangian map. The number of operations can thus be reduced from a 
brute-force estimate O^N"^) to 0(A^ln2 A^). 

In the simulations it is of course necessary to introduce both large-scale 
and small-scale cutoffs. The small-scale cutoff must be chosen sufficiently 
small compared to the integral scale to leave room for a k'"^ inertial range. 
Since the integral scale grows with the time, this constraint is not very strin- 
gent. In practice, we take the initial integral scale about a factor ten larger 
than the mesh. 

For the large-scale cutoff, we assume periodic boundary conditions which 
facilitate Fourier transformations and make no significant difference as long 
as we analyse scales much smaller than the spatial periodicity. 

In our simulations, the maximum wavenumber is taken to be 1/2, so that 
wavevectors run from ±1/A^ to ±1/2, where A^ is the resolution, that is 
the number of grid points in physical space. In all our simulations we take 

= 2^^ ^ 2 X 10^. The energy spectrum for the initial Gaussian velocity is 

E,{k) = al\k\-Q~^^ , A; = 0,±^,±|,...,±^, (122) 

with fco ~ 0.07 and a„ chosen such that {"ipoix)) = 1. The initial integral 
scale and nonlinear time, as defined in (|14D and (^4]) , respectively, are then 
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given by 

1/2 



ko 



r 



n-l 
2 



n.+l 
2 



tni = T2 TTk^TY- (123) 



^0 rf^ 



With the chosen value of ko, the integral scale is around ten and the nonlinear 
time around one hundred. 

The initial potential is generated by fast Fourier transform (FFT) from 
its Fourier coefficients ipo{k) which are complex Gaussian variables of vari- 
ance {\ipo{k)\'^) = Eo{k)/{Nk'^). After the solutions at suitable output times 
have been constructed with the procedure described above, they are Fourier 
transformed and their energy spectra are calculated by averaging typically 
over 10^ realizations.El 

The precise value of the spectral exponent 1 < n < 2 in our runs depends, 
of course, on the particular tradeoffs desired: the closer n is to 1, the more 
conspicuous is the persistence of the initial range; and the closer it is to 
2, the more conspicuous is the (nonlinearly generated) k'^ range. A natural 
unit for the output times is t* = N"^ . Indeed, we know from ( ^Tf ) that the 



integral scale grows roughly as t^^^. Thus t* is roughly the time needed for 
the integral scale to grow from an initial value 0(1) to a value comparable 
to the largest scale available in the simulation, which is 0{N). In practice, 
we take output times 

^i = 10'-^A^^ i = l,2,...,8. (124) 

We describe now the results of our numerical experiments. Fig. ^ shows 
the evolution of the spectrum for n = 1.7 and is mostly meant to bring out 
qualitative aspects. At large wavenumbers a inertial range develops very 
quickly.§ At short times the spectrum preserves its initial value at all but 
the largest wavenumbers, in accordance with the permanence of large eddies. 
The wavenumber where the spectrum achieves its maximum, which we use 
as an operational measure of the inverse integral scale at time moves to 
smaller and smaller wavenumbers. At large times the spectrum appears to 
change in a self-similar way. 



^^This Monte Carlo procedure is easily parallelized in the implementation on an IBM 
SP2 machine. 

^°The time scale for shocks to develop is t^i- 

^^It is actually almost exactly equal to this (see Gurbatov, Malakhov & Saichev 1991, 
Fig. 5.13). 
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wavenumber 



Figure 5: Evolution of the energy spectrum with an initial spectrum (dashed 
line) proportional to k"" {n = 1.7) at small wavenumbers k. Resolution 
= 2^^. Spectra averaged over 10^ realizations. The labels correspond to 
output times U = 10^"^ A^^ = 1, 2, ... 8). 

We turn to a more quantitative description. Fig. ^ shows the energy 
spectrum at the third output time ts = 10~^A^^. We also plot the initial 
spectrum with its \k\"' tail which should be preserved at small wavenumbers 
(leading term of outer expansion) and the leading term of the inner expan- 
sion, obtained by Fourier transformation of (0), which has a fc^ tail at small 
wavenumbers. The intersection of the latter two is the switching wavenumber 
ks{t) as defined by ( [L19| ). It is seen that, except in the immediate neighbor- 
hood of the switching wavenumber, the leading terms of the inner and outer 
expansions are in very good agreement with our simulations. 

Fig. 1^, communicated by A. NouUez, gives additional evidence for the 
permanence of large eddies. 
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Figure 6: Energy spectrum at t = = 10"'^ A^^. Same conditions as Fig. ^ 
Initial spectrum : dashed line. Leading-order asymptotic theory in inner re- 
gion : dashed dotted line. The vertical line shows the switching wavenumber. 



We examine now the temporal evolution of various quantities for which 
we have theoretical predictions. Fig. ^ shows the evolution of the integral 
scale compared to the theoretical prediction(^). The agreement is excellent 
provided that the logarithmic correction is taken into account. By ( ^8]) and 
( |29D this also implies the validity of Kida's log-corrected law for the energy 
decay. Fig. |^ shows the corresponding information for the coefficient A[t). 

Fig. |TD| shows the evolution of the switching wavenumber ks{t), compared 
with the theory (|119|) . Finally, Fig. |Tl| shows the evolution of the switching 
wavenumber, but now for a spectral index n = 1.5. The agreement with the 
leading-order theory is slightly better than for n = 1.7. 
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Figure 7: Evolution of the energy spectrum compensated by a A;~" factor 
in order to reveal the region of "permanence of large eddies" (PLE) as a 
plateau with unit height. Spectral exponent n = 1.5. Resolution = 2^°. 
Spectra averaged over 30,000 realizations. The unit of time is different from 
the other figures, the spatial period being here unity. Observe the progressive 
shrinking of the PLE zone. {Courtesy A. Noullez.) 

9 Concluding remarks 

The main conclusions have already been presented in Section |^ (see also the 
Abstract). Here, we make additional comments. 

We have given an explicit formula only for the leading-order term of the 
law of energy decay. The higher-order terms have just been bounded : they 
appear to be smaller (in relative terms) by negative powers of the time t. It 
would be useful to have explicit expressions for the subdominant terms. 

Our derivation of the law of energy decay uses an asymptotic expression 
for the probability of non-crossing of the initial potential and suitable parabo- 



43 




time 



Figure 8: Evolution of the computed integral scale L{t) (asterisks) compared 
to the theoretical leading-order prediction (solid line) and the same without 
the logarithmic correction (dashed line). Same conditions as Fig. ^. 

las. Kida's (1979) original derivation obtained this from the distribution of 
shocks, which he related to the distribution of (high local) maxima in the 
initial potential. Actually, it should be possible to reformulate our theory 
in terms of maxima rather than of upcrossings. Roughly, we expect that 
high-level upcrossings should occur very close to high-level maxima. This 
statement cannot be exactly true since two successive up- and downcrossings 
of the (initial) potential with a parabola need not be separated by a maximum 
of the potential. However, crossings of a parabola having a high-lying apex 
are very likely to happen in the fiat region near the apex of the parabola 
where up- and downcrossings are separated by a maximum. Quantifying 
such remarks would be of special interest for extension of the theory to the 
multidimensional Burgers equation. Indeed, in > 1 space dimensions, the 
equivalent of the crossings are (rf — l)-dimensional manifolds for which the ex- 
tension of the method presented here is far from straightforward. But (local) 
maxima are still point processes and can be handled via the (i-dimensional 
version of the Balian-Schaeffer formula (Balian & Schaeffer 1989). 
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Figure 9: Evolution of the coefficient A{t) in front of the k"^ term at inter- 
mediate wavenumbers (asterisks), compared to the leading-order theoretical 
prediction (|117|) (solid line) and the same without its logarithmic correction 
(dashed line). Same conditions as Fig. |[ 

There are a number of other interesting extensions of the present work, 
such as the effect of viscosity and of non-Gaussian initial potentials. Instances 
of these have already been studied by Funaki, Surgailis & Woyczynski (1995), 
Surgailis (1996), Hu & Woyczynski (1996) and Surgailis (1997). 

Finally, a completely open question is the extension of the present ideas 
to three-dimensional incompressible turbulence. Does the mechanism which 
prevents globally self-similar decay for 1 < n < 2 have a counterpart for the 
Navier-Stokes problem? 
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Figure 10: Evolution of the switching wavenumber ks{t), compared to the 
leading-order theoretical prediction ( |119D (solid line) and the same without 
its logarithmic correction (dashed line). Same conditions as Fig. |^. 
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Figure 11: Evolution of the switching wavenumber ks{t) for n = 1.5. Other- 



wise, as in Fig. 10 



Appendix 

A The Balian— Schaeffer formula for "parti- 
cle" distributions 

The original work of Balian & Schaeffer (1989) was concerned with the fol- 
lowing type of problem. One has a distribution of "particles" (galaxies, in 
the applications considered by them) in the three-dimensional space which 
is homogeneous, i.e. translation-invariant and one wants to determine the 
probability that there is no particle at all within some finite domain. For the 
application to the Burgers equation, the role of the particles will be played 
by the upcrossings between the initial potential and a given parabola, which 
constitute a point process in the sense of Leadbetter, Lindgren & Rootzen 
(1983). Although the potential is a homogeneous random function, the pres- 
ence of the parabola makes our problem inhomogeneous. Furthermore, our 
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domain is the whole real line and is thus not bounded.c2l Nevertheless, our 
approach will remain rather close in spirit to the original Balian-Schaeffer 
approach. In the remainder of this appendix we shall use the word "particle" 
rather than "intersection", since the formalism has obviously applications 
beyond the Burgers equation. 

Let there be given a random set of particles on the real line. The total 
number of particles is assumed to be random but almost surely finite. We 
are actually interested in finding the probability of having no particle at all. 
We shall characterize this random set by the so-called n-particle densities. 
Fi{yi) dyi is the probability of having a particle within the interval [yi, yi + 
dyi]. Similarly, Fn{yi, y2, ■ ■ ■ , y-n) dyidy2 . . . dyn is the probability of having 
(at least) n different particles within, respectively, the intervals [yi,yi + dyi], 
[z/2,1/2 + dy2], ■ ■ - and, finally, [yn,yn + dUn]- Note that when the particles 
are independent (an assumption not made here), F„ is just a product of Fi 
factors. Note also that the F„s are not normalized. For example / Fi{yi)dyi 
is the mean number of particles. ^1 Balian and Schaeffer's work leads to 
expressions for the probability of having exactly n particles in the domain, 
in terms of the F„s. We shall here concentrate on the evaluation oi po, the 
probability of having no particle at all. 

We begin by introducing the exclusive n-particle densities corresponding 
to the case where there are exactly n particles. Pn{yi,y2, ■ ■ ■ , Un) dyidy2 ■ ■ ■ dyn 
is the probability of having a particle within the interval yi + dyi\, another 
particle within the interval [1/2,1/2 + '^1/2], • • • , a particle within the interval 
[l/n,l/n + c^l/n] and no other particle. Observe that 

j Pniyi,y2,---,yn)dyidy2...dyn = n\pn. (A.l) 

The factor n\ is present because the P„s are here defined for unlabeled par- 
ticles. We shall use ( |A.1| ) to calculate p„ for n = 1, 2, .. and then obtain po 



from the obvious relation 

Po = l-pi-p2- ...Pn- (A. 2) 



^^The formalism developed hereafter applies identically to bounded and unbounded 
domains. The integral sign will always denote a simple or multiple integration over the 
whole domain. 

^■^This follows from the observation that the mean number of particles in the interval 
[yijj/i + dyi] differs by 0{dyi) from the probability of having one particle within the 
interval [yi,yi + dyi]. 
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expressing that the total probabihty is unity. 

We show now that the p„s may be expressed in terms of the -F„s. First 
we relate the F„s to the P^s. To have a particle in [yi, yi + dyi] we must have 
one of the following mutually exclusive situations : (i) a total of exactly one 
particle in [yi,yi + dyi], (ii) a total of exactly two particles, one of them in 
[2/1 5 yi + dyi] and the other one at any place, .... Hence, 

= Pi{yi) + J P2{yi,y2)dy2 

+^ / ?/2, ys) dy2dy3 + .... (A.3) 

The presence of the factorials in the denominators comes again from the 
unlabeled character of the particles. Similarly, we have 

Fn{yi, ■■■yn) =Y.n / Pn+liVl, • • • 2/n, ?/n+l, • • • , yn+l)dyn+l ■ ■ ■ dyn+l- (A.4) 

Next we introduce the quantities 

fn = J Fniyu 2/2, • • • , yn)dyidy2 . . . dyn- (A.5) 
Integrating (|A.3| ) over dyi, we obtain 

h=Pl + 2\p2 + -P3 + ... -, —Pm + .... (A.6) 

2! (m — 1)! 

More generally, integrating ( |A.4| ) over dyi, . . . , dyn, we have 

fn = n\pn +{n+ l)\pn+l + • • • 7 '^,Pm + .... (A. 7) 

[m — n)\ 

Multiplying ( [A.7| ) by (— )"/?t,! and summing over n, we obtain 

Po = E -^fn- (A.8) 



n=0 



Let us briefly digress and consider the case of independent particles, when 

Fn{yi, ...,yn) = Fi{yi)Fi{y2) ■ ■ ■ Fi (?/„), and thus /„ = . From (^) we 
then obtain 

Po = exp (-/i) . (A.9) 
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Since fi = J Fi{yi)dyi is the mean number of particles, ( |A.9| ) may be viewed 
as just a consequence of the Poisson distribution. It is referred to in this 
paper as the "Poisson approximation" . 

The general case, of correlated particles, is handled by introducing corre- 
lation functions through a cluster expansion : 



^2(2/1,1/2) = C'i(?/i)Ci(2/2) +(^2(2/1,2/2) 



(A.IO) 
(A.ll) 



Fn{yi,y2,- ■ ■) = ■■■Ciiyri,---,yn) 

all partitions 

We also need the cumulants 

Cn = y Cniyi, 2/2, ■ • • , yn)dyidy2 . . . dyn- 



(A.12) 



(A.13) 



Note that the c„s for n > 1 vanish in the case of independent particles. 
The name "cumulant" for the c„s is justified by the observation that the c„s 
and the /^s defined in ( |A.5|) are related like moments and cumulants: if we 
introduce the generating function 



00 



we easily obtain from (|A.1CI| )- (|A.12|) 

/ 00 

n{z) = exp ^ 



(A.14) 



(A.15) 



Expansion of the r.h.s. of ( |A.15 ) in powers of z and identification with the 
expansion ( [A.14|) allows the successive determination of the set of /„s in 
terms of the set of c„s and, conversely, identification of the expansion of the 
logarithm of the r.h.s. of ( |A.14|) with the logarithm of ( |A.15|) gives the c„s 
in terms of the /„s. In particular, we have 

Cl = fi = J Fi{y,)dy^, (A.16) 

C2 = /2 - fi = J (^2(2/1, 2/2) - Fi(2/i)i^i(2/2)) dy,dy2. (A.17) 
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Comparing ( |A.8| ) and ( |A.14| ), we have po = Q{—1). Then, using ( |A.15| ), 
we obtain the required expression for the probabihty of having zero particles : 



Of course, if we truncate the series to its first term, we recover the Poisson 
approximation. 
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FIGURE CAPTIONS 



Figure 1: Sketch of the energy spectrum at long times when 1 < n < 2. 

Figure 2: Intersections of the initial potential ipo^y) (curve Gi) with the parabola 
G2. 

Figure 3: Intersections of the initial potential ipQ with two parabolas. 

Figure 4: The different regions of integration appearing in ( |106D . Near the diag- 
onal Hi = y2 we have regions I and I' of ( |107|) . Away from the diagonal 
we have regions II and 11'. 

Figure 5: Evolution of the energy spectrum with an initial spectrum (dashed line) 
proportional to /c" [n = 1.7) at small wavenumbers k. Resolution = 
2^^. Spectra averaged over 10^ realizations. The labels correspond to 
output times U = 10''^ N^, (i = 1, 2, ... 8). 

Figure 6: Energy spectrum at t = ^3 = lO^'^A^^. Same conditions as Fig. ^ 
Initial spectrum : dashed line. Leading-order asymptotic theory in inner 
region : dotted line. The vertical line shows the switching wavenumber. 

Figure 7: Evolution of the energy spectrum compensated by a /;;"" factor in order 
to reveal the region of "permanence of large eddies" (PLE) as a plateau 
with unit height. Spectral exponent n = 1.5. Resolution A^ = 2^°. 
Spectra averaged over 30,000 realizations. The unit of time is different 
from the other figures, the spatial period being here unity. Observe the 
progressive shrinking of the PLE zone. [Courtesy A. Noullez.) 

Figure 8: Evolution of the computed integral scale L{t) (asterisks) compared to 
the theoretical leading-order prediction (solid line) and the same with- 
out the logarithmic correction (dashed line). Same conditions as Fig. ^ 

Figure 9: Evolution of the coefficient A{t) in front of the k"^ term at intermediate 
wavenumbers (asterisks), compared to the leading-order theoretical pre- 
diction ( p.l7| ) (solid line) and the same without its logarithmic correction 
(dashed line). Same conditions as Fig. ||. 

Figure 10: Evolution of the switching wavenumber /cs(t), compared to the leading- 
order theoretical prediction ( |119| ) (solid line) and the same without its 
logarithmic correction (dashed line). Same conditions as Fig. ^ 

Figure 11: Evolution of the switching wavenumber ks{t) for n = 1.5. Otherwise, 
everything as in Fig. 0. 
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